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ABSTRACT 

Context. The tidal evolution of planets orbiting brown dwarfs (BDs) presents an interesting case study because BDs' terrestrial planet 
forming region is located extremely close-in. In fact, the habitable zones of BDs range from roughly 0.001 to 0.03 AU and for the 
lowest-mass BDs are located interior to the Roche limit. 

Aims. In contrast with stars, BDs spin up as they age. Thus, the corotation distance moves inward. This has important implications 
for the tidal evolution of planets around BDs. 

Methods. We used a standard equilibrium tidal model to compute the orbital evolution of a large ensemble of planet-BD systems. We 
tested the effect of numerous parameters such as the initial semi-major axis and eccentricity, the rotation period of the BD, the masses 
of both the BD and planet, and the tidal dissipation factors. 

Results. We find that all planets that form at or beyond the corotation distance and with initial eccentricities smaller than ~ 0.1 are 
repelled from the BD. Some planets initially interior to corotation can survive if their inward tidal evolution is slower than the BD's 
spin evolution, but most initially close-in planets fall onto the BD. 

Conclusions. We find that the most important parameter for the tidal evolution is the initial orbital distance with respect to the 
corotation distance. Some planets can survive in the habitable zone for Gyr timescales, although in many cases the habitable zone 
moves inward past the planet's orbit in just tens to hundreds of Myr. Surviving planets can have orbital periods of less than 10 days 
(as small as 10 hrs), so they could be observable by transit. 

Key words. Brown dwarfs - Stars: rotation - Planets and satellites: dynamical evolution and stability - Planet- star interactions - 
Astrobiology 



1. Introduction 

Tidal effects around brown dwarfs (BDs) have received lit- 
tle attention but they present an interesting case study be- 
cause the planets are thought to form very close to their star 
in a region where tides are very strong and lead to signifi- 
cant orbital changes and extreme internal planetary heating (see 
Papaloizou & Terquem, 2010; Barnes et al., 2010). Moreover, as 
the luminosity of brown dwarfs decreases dramatically on Gyr 
timescales (see Tout et al., 1999), the habitable zone moves in- 
ward in time and meets the formation region of planets, which 
naturally raises the question of habitability of these planets. 

The parameter that governs the direction of tidal evolution 
for a planet orbiting a star or BD (or a satellite orbiting a planet) 
is the initial semi-major axis with respect to the corotation ra- 
dius, the orbital radius where the orbital period matches the cen- 
tral body's spin period. For a planet interior to the corotation 
radius, the planet's mean motion is slower than the primary's ro- 
tation so the tidal bulge raised by the planet on the primary lags 
behind the position of the planet. The planet feels a drag force 
that slows it down and causes its orbital radius to shrink, in some 
cases leading to an eventual merger with the primary. However, 
for a planet exterior to the corotation radius, tidal forces push the 
planet outward because the tidal bulge on the star is in advance 
with respect to the position of the planet. 

Main sequence stars spin down in time due to angular mo- 
mentum loss driven by stellar winds (Skumanich, 1972; Matt 



et al., 2010). In contrast, BDs spin up with time (see Herbst et al., 
2007; Scholz & Eisloffel, 2005). This key difference, which is 
due to the low efficiency of the magnetic braking in the substel- 
lar domain (see Section 2.3.3), is of great importance for the tidal 
evolution of planets orbiting BDs because the primary's spin rate 
sets the position of the corotation radius. For stars, the corotation 
distance expands as the star spins down, so close-in planets spi- 
ral inward and, in time, more distant planets also start to spiral 
inward. However, the corotation distance around BDs shrinks in 
time such that some planets may spiral inward at early times but 
later, when the corotation radius moves interior to their orbits, 
change direction and spiral outward, away from the BD. 

The habitability of planets around BDs also differs from 
planets around main sequence stars. Given the inability of BDs 
to fuse H into He, BDs simply cool in time and their habitable 
zones move inward. For a planet on a fixed orbit, the habitable 
zone will thus sweep past it in a finite time (Andreeshchev & 
Scalo, 2004). For a planet being pushed outward by tides, the 
"habitable interval" will be even shorter. The time a planet will 
spend in the habitable zone depends on many parameters, the 
main one being the mass of the star. 

Furthermore, BDs small sizes and masses also make them 
excellent targets for radial velocity and transit searches (see 
Blake et al., 2008), although significant observational constraints 
exist (e.g., the difficulty of measuring precise radial velocities in 
the infrared (see Blake et al., 2010)). BDs are thought to be more 
numerous than low mass stars, and recently five new BDs have 
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been discovered at ~ 10 pc (Burgasser et al., 2011) and two at 
~ 5 pc (Scholz et al., 2011) with WISE increasing the number 
of targets for observation of planetary transit. 

There are several reasons to think that planets should form 
around BDs. The disk fraction of young BDs is comparable to 
that for T Tauri stars (see Jayawardhana et al., 2003; Luhman 
et al., 2005), and grain growth and settling have been observed 
in disks around BDs (Apai et al., 2005). Protoplanetary disk 
masses are thought to correlate roughly linearly with the stel- 
lar mass (Scholz et al., 2006; Andrews & Williams, 2007), such 
that BD disks might be systematically lower-mass than those 
around main sequence stars. However, recent sub-mm observa- 
tions suggest that there exists a correlation between protoplane- 
tary disk mass and radial extent such that lower-mass disks are 
more compact (Andrews et al., 2010), and disks around low- 
mass stars may actually last longer than around more massive 
stars (Pascucci et al., 2008). In addition, Kepler observations 
have found that low-mass stars actually harbor close-in super 
Earth- sized planets at a higher frequency than Sun-like stars 
(Howard et al., 201 1). Thus, we expect that there probably exists 
a sizable niche in the parameter space of BD and disk properties 
(BD mass, disk mass, surface density distribution, and lifetime) 
that allows for the formation of sizeable (~ 0.1 Earth masses 
or larger) terrestrial planets near the habitable zone of BDs (see 
Raymond et al, 2007; Payne & Lodato, 2007). 

Here we use a standard equilibrium model to study the tidal 
evolution of planets orbiting BDs. Our paper is structured as fol- 
lows. The tidal and BD evolutionary models are discussed in 
Section 2. The results of the tidal evolution of planets around 
BDs and the influence of several parameters are presented in 
Section 3. Finally, in Section 4 we discuss tidal evolution with 
regards to the habitable zone around BDs and the importance of 
tidal evolution for observations. 



and depends on the mass of the planet M p , its dissipation cr p and 
of the mass of the BD M*. Q p is the planet rotation frequency, 
and n is the mean orbital angular frequency. The stellar param- 
eters are obtained by switching the p and * indices. Nal(e) and 
Na2(e) are eccentricity-dependent factors, which are valid even 
for very high eccentricity (Leconte et al., 2010) : 



Nal(e) = 
Na2(e) = 



1 + 31/2e 2 + 255/8e 4 + 185/16<* 6 + 85/64<* 8 

(1 - e 2 ) 1 ^ 2 
1 + 15/2e 2 + 45/8e 4 + 5/16<? 6 
(1 - e 2 ) 6 ' 



Similarly, the evolution of the eccentricity of the planet is 
given by : 
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Finally, the rotational state of each object is given by : 
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2. Model description 

We have developed a model to study the orbital evolution of 
planets around BDs by solving the tidal equations for arbitrary 
eccentricity and also taking into account the observed spin evo- 
lution of BDs. 



and 
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2.1. Tidal model 

The tidal model that we used is a re-derivation of the equilib- 
rium tide model of Hut (1981) as in Eggleton et al. (1998). 
We consider both the tide raised by the BD on the planet and 
by the planet on the BD. We use the constant time lag model 
(Leconte et al., 2010) and use the internal dissipation constant cr 
which was calibrated for giant exoplanets and their host stars by 
Hansen (2010). 

The secular tidal evolution of the semi-major axis a is given 
by (Hansen, 2010): 



1 da 
a dt 



1 r l 
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(1) 



where the dissipation timescale T p is defined as 



1 



1 



9M*(M v + M*)RM<r p 



(2) 



Here, y t = is the ratio of orbital angular momentum h to 
spin angular momentum and rg2* is the square of the parameter 
rg* (which is the radius of gyration of Hut (1981)) which is de- 
fined as : / = M*(rg*R*) 2 , where / is the moment of inertia of 
the BD, and 



Nol(e) = 



No2(e) = 



l + 15/2e 2 +45/8e 4 + 5/16e 6 
(1 - e 2 ) 13 / 2 ' 
1 + 3e 2 + 3/8e 4 
(1 - e 2 ) 5 ' 



The last two terms of equation (5) translate the conservation 
of the moment of inertia. As will be discussed in the next sub- 
section, in the case of BDs the radius decreases with time. 

The integration of these equations was performed using a 
fourth order Runge-Kutta integrator with an adaptive timestep 
routine (Press et al., 1992). The precision of the calculations was 
chosen such that the final semi-major axis of each integrated sys- 
tem was robust to numerical error at a level of at most one part 
in 10 3 . 
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2.2. Planets model 

Given the relatively small mass of protoplanetary disks around 
BDs, we assume that most planets will be terrestrial planets 
rather than gas giants. Although we do not exclude the possi- 
ble existence of gas giants around BDs (see Section 3.2.5), we 
focus on terrestrial planets whose internal dissipation factor is 
the same as Earth's but scaled with the radius of the planet. 

Neron de Surgy & Laskar (1997) inferred the quantity 
k 2 ,®At (B = 213 s from the DE245 data for Earth. k 2 , p is the Earth's 
potential Love number of degree 2, which is a parameter depend- 
ing on the moment of inertia of the body. It tells us how the body 
responds to compression (k 2 = 3/2 means that the body is an 
incompressible ideal fluid planet). AT p is the constant time-lag. 
Hansen's cr p and the quantity k 2 , p AT v are related through : 



k 2 , p AT p = 



3 gg> 
2 Q ' 



(6) 



where R p is the planetary radius and Q is the gravitational 
constant. 

The planet's compositions are assumed to be an Earth-like 
mixture of rock and iron, following the mass-radius relation of 
Fortney et al. (2007): 



RJR @ = (0.0592 rmf + 0975) (log^f 

+ (0.2337 rmf + 0.4938) log— p - 
+ (0.3102 rm/ + 0.7932), 



(7) 



where rmf is the rock mass fraction. For Earth, rmf ~ 0.7. 

Concerning the evolution of the planet it is interesting to note 
that Equation (4) has a much shorter timescale than the other 
differential equations. Very quickly in its evolution the planet 
reaches pseudo-synchronization, which means that its rotation 
tends to be synchronized with the orbital angular velocity at peri- 
astron (Hut, 1981). Indeed, the tidal interaction is strongest near 
periastron where the planet is closest to the BD. Having verified 
this behavior in preliminary simulations, the rotation period is 
fixed to the pseudo- synchronization value at every calculation 
timestep. 

In this study, we chose not to treat the evolution of the obliq- 
uity of the bodies for simplicity. Heller et al. (201 1) showed that 
the timescale of evolution of obliquity - which they call "tilt ero- 
sion" time - is shorter for low mass stars and for close-in planets. 
BDs have low masses and the planets we consider in this study 
are very close-in, so we can assume that the obliquity is negligi- 
ble. 

2.3. Brown Dwarf Evolution 

BDs are self-gravitating objects that formed through gravita- 
tional collapse of a dense molecular cloud but that are not hot 
enough to ignite proton fusion (REF BD review). As the hydro- 
gen burning temperature is around 3 x 10 6 K, only objects more 
massive than M H bmm ~ 75 M JUP are able to start the PPI nu- 
clear reaction chain, providing us with an upper bound for the 
BD mass domain (Chabrier & Baraffe, 1997, 2000). If the lower 
bound is more uncertain, both observational (Caballero et al., 
2007) and analytical (Padoan & Nordlund, 2004; Hennebelle & 
Chabrier, 2008) arguments suggest that the same formation pro- 
cess responsible for star formation can continuously produce ob- 
jects down to a few Jupiter masses. Mini brown dwarfs and giant 



planets (formed in a protoplanetary disk) thus overlap in mass, 
stressing the need for identification criteria enabling the distinc- 
tion between these two types of astrophysical bodies (Leconte 
et al., 2009; Spiegel et al., 2011; Leconte et al., 2011a). 

2.3.1 . Internal structure and luminosity 

To calculate the influence of the dissipation into the BD on the 
tidal evolution of the system, one needs to know the internal 
structure of the object and how it responds to a tidal perturbation, 
namely, the BD radius, moment of inertia (the gyration radius), 
tidal response (the love number, k 2 ) and tidal dissipation (<x). 
The BD radius, moment of inertia and tidal response are pro- 
vided by the grids of BDs and giant planets evolutionary models 
of Leconte et al. (201 lb) 1 . While the main physics inputs (equa- 
tions of state, internal composition, atmosphere models, bound- 
ary conditions) used in these calculations have been described 
in detail in previous papers devoted to the evolution of substellar 
objects (Chabrier & Baraffe, 1997; Baraffe et al, 2003a; Leconte 
et al., 2011b) details the calculation of the dynamical properties 
of gaseous objects in hydrostatic equilibrium (inertia and love 
number). 

Because most of the BD's interior remains fully convective 
for most of its lifetime, there is no drastic change in the dimen- 
sionless moment of inertia of the object with time, as is the case 
in young solar type stars when the radiative core develops. The 
main variation of the moment of inertia of the BD is due to its 
monotonous contraction caused by the radiative energy losses. 
However, if the object is more massive than ~ 0.0125 M© the 
temperature reached in the interior due to the early contraction 
is sufficient to ignite Deuterium burning (this limit is ~ 0.06 M© 
for Lithium burning). As this additional energy source is tem- 
porarily able to nearly balance the radiative losses, the contrac- 
tion is slowed down and the radius, the luminosity, and the mo- 
ment of inertia reach a plateau. Because the primordial abun- 
dance of these compounds is small, this phase lasts only up to a 
few million years for massive brown dwarfs, and 100 Myr near 
the deuterium burning limit. 

To quantify the extent of the habitable zone we use a simple 
criterion based on the flux at the planet's substellar point. At a 
given semi major axis and eccentricity, this mean bolometric flux 
is 



eff 



(8) 



where <tsb is the Stefan Boltzmann constant, and R B d, and T e s 
are the radius and effective temperature of the BD. 

2.3.2. BD dissipation 

Hansen (2010) makes the assumption that all substellar objects 
have the same internal dissipation because BDs and gas giants 
have similar internal structure and dissipative response. Note, 
however, that if the presence of a central dense (not necessarily 
solid) core plays a major role in the process responsible for the 
tidal dissipation in giant planets (Goodman & Lackner, 2009), 
this assumption may not hold anymore, as BDs, which form like 
stars are not expected to possess a metal enriched core. In our 
notation, the internal dissipation factor ct B d is thus the same as 
cr p from Hansen (2010). 



1 Electronic versions of the model grids are avail- 
able at http://perso.ens-lyon.fr/jeremy.leconte/JLSite/JLsite 
/Exoplanets Simulations, html 
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The dissipation factor is poorly constrained but Hansen 
(2010) gives estimates for stars : <x* = 6.4 x 10" 59 g" 1 cm" 2 s" 1 o = * 
where cr* = 7.8 x 10" 8 and for gas giants : cr p = 5.9 x 
10" 54 g" 1 cm" 2 s" 1 o = /? where ^ = 3.4 x 10" 7 . Thus for stars, the 

and for gas 



Radius and rotation period of BDs of different masses 



-1 -2 - 

cm z s 



-1 -2 - 

cm z s 



dissipation factor is <x* = 4.992 x 10~ 76 
giants or BDs it is cr p = 2.006 x 10" 60 g" 

So the value of the dissipation factor for the BDs thereafter 
referred as <tbd is : 



cr BD = 2.006 x 10" 53 kg- A m- z s" 



-l_-2^-l 



(9) 



2.3.3. Rotational angular velocity 



Brown dwarfs are very rapid rotators, which is surprising at 
first glance. The rapid rotation coupled to the convective interior 
creates a strong magnetic field which should entail an impor- 
tant magnetic braking through ejection of a magnetized wind. 
However, as shown by Mohanty et al. (2002), in BDs' cool at- 
mospheres the degree of ionization is very low. A combination 
of very low ionization fraction and high densities results in very 
large resistivities and thus efficient magnetic field diffusion. The 
magnetic braking is thus very inefficient in the substellar domain 
and cannot counteract the early spin up caused by the contrac- 
tion. Thus, throughout this study, magnetic braking is neglected. 

Information about the rotation periods of BDs comes from 
both statistical studies of BDs in open clusters (see Herbst et al., 
2007, and references therein), for which ages are known, and 
specific studies of particular BDs, which can be studied in detail 
but for which age determinations are often challenging. 

BDs have observed orbital periods as short as a few hours 
(Herbst et al, 2007). For example, the nearby BD DENIS-P 
J0041 353-5621 12 has a rotation period of 2.8 hr, an age between 
10 and 15 Myr and a mass below 0.062 M (Koen, 2010). The 
rotation period depends on the age of the BD and on its mass 
(Herbst et al., 2007). For a given age, low-mass BDs rotate faster 
than high-mass BDs, and all BDs spin slower at younger ages. 
This behavior can be explained by the simple conservation of 
moment of inertia during BD contraction. 

The assumed evolution of the spin period of BDs adopted 
in our simulations is shown in Fig. 1. After 1 million years, 
BD rotation periods are distributed between 8 hr and 70 hr with 
longer initial periods for more massive BDs (the 70 hr period 
corresponds to the highest mass BD of 0.08 M ). The fact that 
higher-mass BDs have slower initial spin rates is an empirically- 
calibrated relation (Herbst et al., 2007), although the data are 
relatively scarce. We chose these initial values to match the or- 
der of magnitudes seen in Herbst et al. (2007) and we ensured 
that the final rotation periods were not too close to the break up 
period of the BD (i.e., the Keplerian period of a massless par- 
ticle at the surface of the BD). Considering the uncertainty on 
the dissipation factor and the strong gravity at the surface, we 
ignore the increase of the mean radius of the BD caused by the 
fast rotation (see Fig. 1 of Leconte et al. (2011b)). 

We note that these initial conditions remain somewhat un- 
certain because the BD radii are uncertain at young ages. 

In our calculation we took into account the contraction of the 
BD, as seen in Equation 5. For convenience, this equation can be 
partially integrated and computed as follows : 



rg2 BD (t ) I R BD (t ) 
rgl BD (t) \ R BD (t) 



f 

J t 



Xexp ftidesdt 




Fig. 1. Spin evolution of BDs assumed in our simulations, calibrated 
to match observations (Herbst et al., 2007). The curves represent the 
rotation period of BDs of different masses (0.01, 0.012, 0.015, 0.02, 
0.03, 0.04, 0.05, 0.06, 0.07, 0.072, 0.075, 0.08 M ) computed from the 
radius evolution and the conservation of angular momentum. The dotted 
lines represent the initial time for the two sets of simulations, thus also 
giving the initial rotation periods for each BD. 



where to is the initial time and f ti des is given by : 
1 d£l BD I 



ftides — 
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~[Nol(e) 
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No2(e)\. 



(10) 



For terrestrial planets, f^ es remains small throughout the 
evolution, and the BD rotation period is mainly determined by 
the conservation of angular momentum and so by the initial ro- 
tation period. 

3. Results 

3.1. An example 

We now present an example of the tidal evolution of planets 
around a BD. We consider a BD of mass 0.04 M©, and a planet 
of 1 M e placed at initial orbital radii between 7 x 10" 4 AU 
and 2.4 x 10" 3 AU. The BD dissipation is the same as used 
by Hansen (2010) as given in equation (9). The planet has the 
same k2, p At p ans cr p as Earth. The initial BD rotation period is 
36 hr and the initial orbital eccentricity is 0.01. The initial time 
is 10 6 yrs, meaning that the BD's structure at the start of the 
integration was taken from a 1 Myr snapshot in the BD mod- 
els described above. This means that we consider that the planet 
was fully-formed at the time of the dissipation of the disk, which 
is consistent with simple estimates of the accretion time around 
very low-mass stars (Raymond et al., 2007). Of course, given 
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that gaseous protoplanetary disks have lifetimes of a few Myr 
(see Haisch et al., 2001; Hillenbrand et al., 2008) and that disks 
around low-mass stars tend to be longer-lived than around Sun- 
like stars (Pascucci et al., 2009), values for "time zero" as long 
as 10 Myr are probably realistic (see Section 3.2.4). 

The range of evolutionary pathways for this example is 
shown in Figure 2. As other parameters are held fixed, the evo- 
lution is determined by the initial position of the planet with re- 
spect to the corotation distance. 

If the planet's mean motion is slower than the BD's rotation 
rate - ie if the planet's initial semi-major axis is larger than the 
corotation radius - then the tidal bulge raised by the planet on the 
BD is in advance with respect to the position of the planet. The 
resulting effect of that misalignment is that the planet is acceler- 
ated and its semi-major axis increases. Eventually, as the radius 
of the BD decreases the tidal dissipation in the BD decreases and 
the planet evolution stops at a given semi-major axis. This final 
position depends on a number of parameters as will be presented 
in the next subsection. 

On the other hand, if the planet's mean motion is faster 
than the BD's rotation rate - ie if the initial semi-major axis is 
smaller than the corotation radius - then the tidal bulge raised 
by the planet on the BD lags behind the position of the planet. 
The planet is therefore slowed down and its semi-major axis de- 
creases. If the planet spirals in very quickly then it will directly 
hit the BD surface, but if its evolution is slower - if it survives 
the first few million years - then it will reach the Roche limit 
and be ripped apart. Such planetary destruction is likely to form 
a debris disk around the BD, which is potentially observable with 
next-generation instruments such as JWST. 

A planet's fate depends on the first few million years of evo- 
lution, and its semi-major axis is stabilized in about a hundred 
million years. 

The initial orbital radius offers a first order approximation 
of planets' tidally-driven orbital evolution, but as the BD's spin 
rate increases in time some intermediate cases can be seen. For 
example, if a planet's initial orbital distance is inside the coro- 
tation radius but sufficiently close to it, the planet may undergo 
a two-phase evolution. First, the planet will begin falling onto 
the BD. However, if the BD's rotation period decreases faster 
than the planet's orbit shrinks, then the corotation distance can 
pass by the planet's orbit. Then, the direction of tidal evolution 
is reversed and the planet migrates outward. 

Figure 3 shows the evolution of the rotation periods of the 
BD and planet for the planets initially close to the corotation dis- 
tance in Fig. 2. As described above, a planet that forms inside the 
corotation radius may survive tidal evolution. It begins to fall to- 
wards the BD and its rotation period decreases. However, the BD 
is contracting and its rotation period also decreases. Depending 
on the rate of decrease of the two quantities the BD rotation pe- 
riod may catch on the planet rotation and reverse the direction of 
the planet's tidal evolution. Once the rotation period of the planet 
becomes longer than the BD's, tidal forces push the planet away. 
This intermediate evolution is very rare because the correspond- 
ing range of initial semi-major axes is very narrow. In the case of 
Fig. 2 for a BD of mass 0.04 M© and a planet of 1 M®, the width 
of this range is only ~ 10" 3 AU (from ~ 7.5 to 8 x 10" 3 AU). 

Fig 3 shows that, for the most extreme cases, the planet's in- 
fluence acts to accelerate the BD just before crossing the corota- 
tion radius. For less extreme cases, the planet does not influence 
the BD rotation. 




time (yrs) 

Fig. 2. Evolution of the semi-major axis of a 1 M e planet orbiting a 
0.04 M Q BD. The solid colored lines are for different initial semi-major 
axis. The dashed dotted line represent the radius of the BD and the 
dotted line the corotation radius. The habitable zone is also plotted using 
two different conditions of habitability (solid line : the planet receives 
an incoming flux of 400 W/m 2 , dashed line : 240 W/m 2 , see section 4.1 
for more details). The bold dashed line represents the Roche limit. 



Rotation period 

60 p i i j i 




10 3 10 4 10 5 10 6 10 7 10 8 10 9 10 10 

time (yrs) 



Fig. 3. Evolution of the rotation period of a 1 M e planet evolving around 
a 0.04 M Q BD. The full lines represent the evolution for different initial 
semi-major axis. The dashed lines represent the rotation period of the 
BD. 



3.2. Influence of parameters 

In this section we investigate the effect of each parameter while 
keeping the others fixed at arbitrarily-chosen default values. 
Table 1 lists these default values. 



*0 


M BD 


M p 


ctbd 




eo 


(Myr) 


(M ) 


(M e ) 


(kg^nrV 1 ) 


(kg^nrV 1 ) 




1 


0.04 


1 


2.006 x 10" 53 


8.577 x 10- 43 


0.01 



Table 1. Table of default values of parameters. 



For every set of parameters several simulations were done 
with different initial semi-major axis. 
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3.2.1 . Influence of initial eccentricity 

The initial orbital eccentricity can play a role in the early stages 
of tidal evolution. For an eccentric orbit, the tides raised by the 
BD in the planet are strong and dissipation in the planet is im- 
portant. This leads to a quick decrease of the eccentricity to min- 
imize the planet's internal stress. Once the eccentricity becomes 
negligible, the dissipation in the planet is also negligible and the 
evolution is governed by the tides created in the BD. 

For a fixed semi-major axis, a higher eccentricity implies a 
faster planetary rotation because the planet synchronizes at pe- 
riastron where the orbital velocity is larger (Hut, 1981). Thus, 
the initial semi-major axis that divides behaviors (in- 

ward migration and crash) and & (Inward migration but survival 
of the planet or outward migration) is farther from the BD for 
larger initial eccentricity. 




Fig. 4. Evolution of semi-major axis, eccentricity and rotation period of 
a 1 M e planet evolving around a 0.04 M© BD. Each solid line represents 
a planet with a different initial eccentricity. The initial semi-major axis 

for each Case IS d-switchi 

which increases with eccentricity. The dashed- 
dotted line represent the corotation radius and the dashed line the ro- 
tation period of the BD. The vertical lines show the time at which the 
eccentricity effects become negligible for each initial eccentricity. 



Figure 4 shows the evolution of a suite of planets start- 
ing at a SW i tc h with initial eccentricities of 0, 0.01, 0.05, 0.1, 
0.15, 0.2, 0.3, 0.5, 0.6 and 0.8. The early evolution depends 
strongly on the initial eccentricity. At early times an eccentric 
planet rotates faster than the BD such that the planet's semima- 
jor axis decreases and its eccentricity is damped. The planet's 
orbit shrinks but, due to the eccentricity-dependence of pseudo- 
synchronization, the rotation slows. Then, when the BD's rota- 
tion period drops below the planet's, the planet is pushed away. 
It is interesting to note in Fig. 4 that, for these cases which all 
started at a SW i tc h, the evolution of the different cases converges 
once the eccentricity drops below ~ 0.08. 

For circular planetary orbits the only tidal bulge is the one 
raised on the BD by the planet - referred to as the BD tide - 




Fig. 5. Evolution of eccentricity and of the ratio of planetary over BD 
tides of a 1 M e planet evolving around a 0.04 M Q BD. The full colored 
lines represent the evolution for different initial eccentricity. For a non 
zero eccentricity, the planetary tide dominates. When the eccentricity 
reaches 0.08, the strength of the BD tide is equal to the strength of the 
planetary tide. The vertical lines show the time at which the BD tide 
starts to win over planetary tide. 



and the evolution is controlled by dissipation within the BD. 
However, for eccentric orbits there is an additional bulge raised 
on the planet, called the planetary tide, and the evolution depends 
on dissipation within both the planet and BD. 

If the initial eccentricity is 0, then equation (1) can be sim- 
plified as : 



1 da l_r _ £Ibd ~i 

a dt T BD L n J' 



(12) 



Q p at 



If the orbit of the planet is circular, the planet is in perfect 
synchronization so Q p = n, and the first term of equation (1) 
is zero. As expected, only the tides raised by the planet in the 
BD remain. If Q p = n > ^ < and the planet will fall 

towards the BD (behavior Ji) . If Q p = n < Q BD , ^ > and the 
planet will be pushed away (behavior &). As discussed above, 
the fact that £Ibd increases with time allows for the existence of 
the intermediate case in which the planet starts falling but is later 
pushed outward when £l BD becomes larger than Q p = r 

If the initial eccentricity is non zero, and for Qbd 
time zero, the difference between the strength of BD and plan- 
etary tides comes from the dissipation timescales. When the ec- 
centricity is large, the planetary tide - i.e., the rate of energy dis- 
sipation within the planet - dominates the evolution because the 
amplitude of the tidal bulge varies over the course of the planet's 
orbit. For a 1 M© planet and a 0.04 M Q BD the ratio of the dis- 
sipation timescales is : T BD /T p ~ 10 at to = 10 6 yr and ~ 10 9 at 
1 Gyr. Rocky planets dissipate more than BDs and the difference 
in time is due to the shrinking radius of the BD, which makes its 
dissipation decrease - see equation (2). So while the eccentricity 
is non zero, the tides raised on the planet by the BD will pre- 
vail. In this case, Nal(e) » Na2(e), so ^ is negative and the 
planet is pulled inward. This is always true because for a planet 
in pseudo-synchronization the planetary tide always acts to de- 
crease the orbital distance (see Leconte et al., 2010). But when 
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the eccentricity is small the BD tide dominates and determines 
the evolution of the planet: if the planet is beyond the corotation 
radius, ^ is positive so the planet is pushed outward, and if it is 
inside, ^ is negative so the planet is pulled inward. 

Figure 5 shows the ratio of the strength of BD to planetary 
tides for the case shown in Fig. 4, where the strength of each 
tide is the absolute value of its contribution to ^. While e > 
0.08, the planetary tide dominates but when e < 0.08 the BD 
tide dominates. This is true for all cases studied and analytically 
we can show that for an eccentricity of this order of magnitude 
BD and planetary tides have the same strength for the orbital 
distance evolution. At early times the planetary tide acts to damp 
the orbital eccentricity and shrink the planet's orbit but when the 
eccentricity drops below ~ 0.08 the BD tide takes over and the 
orbit expands. At this point the planet has a constant deformation 
and the ratio of the two tides is very high (> 10 4 ) because the 
planetary tide is so small. As the planet is pushed outward the 
BD contracts and the BD tide also decreases to zero. 



3.2.2. Influence of the BD mass 



Delta a/a 



10 




0.01 0.02 0.03 0.04 0.05 0.06 0.070.08 

Stellar mass (M sun ) 



Fig. 6. Map of Aa/a for a range of initial semi-major axis and different 
BD masses. Time zero is 10 6 yrs. The deep purple zone is the region for 
which the planets don't survive and the white zone is where the planets 
approach the BD without falling on it. The dashed dotted line is the 
corotation distance at the beginning of integration and the dashed line 
is the Roche limit. 



The BD mass is important for tidal evolution because the more 
massive the BD the stronger the tides raised on the planet. The 
BD mass affects its response to the planet's gravity, i.e., the tides 
that the planet raises on the BD (referred to as the BD tides). 

The dissipation timescale for the BD is simply calculated by 
switching the p and * indices in Eqn (2) - and replacing * by BD 
- to yield: 

_ 1 M BD a 8 1 
TBD ~ 9M p (M p + M BD )R)^^ (13) 



3.2.3. Influence of the initial BD spin rate 

Given that there remains considerable uncertainty in the spin 
evolution of BDs, we performed a suite of tidal evolution cal- 
culations assuming that all BDs start with the same initial rota- 
tion period of 36 hours. This exercise isolates the dependence 
on the BD mass, although in a regime that is probably not realis- 
tic given current observations (Herbst et al., 2007). In addition, 
we did not take care in these integrations to avoid very fast final 
rotation periods that might be non physical because they may 
approach breakup velocity. 



The BD mass has an indirect influence on T BD via the BD ra- 
dius R B d- More massive BDs have larger radii - for ages younger 
than 10 8 yrs - and thus heavier BDs have a shorter dissipation 
timescales, which means that tides are stronger. 

The BD mass also has an influence on its rotation, as seen in 
Fig. 1. The initial conditions taken in this study are that heavier 
BDs rotate slower, which is consistent with observations (Herbst 
et al., 2007). The corotation radius is therefore farther from heav- 
ier BDs, and they rotate slower so the tides will be less efficient 
than for lighter BDs, which have smaller corotation radius and 
higher rotation velocities. 

This dependance of the BD mass can be seen in Fig. 6 and in 
Fig. 8 in the next subsection 3.2.4. 

Here, we define the quantity Aa/ao, the fractional change in 
semi-major axis over the planet's evolution, as Aa/ao = (cif - 
0o)M)> where ao and a$ are the semi-major axis at t - to and at 
t = tf = 10 10 yrs, respectively. When this quantity is positive, 
the planet is pushed away and when it is negative it comes closer 
to the BD. When a planet falls onto the BD, Aa/ao = -1. 

In Figure 6, the dark zone is the region for which the planet 
falls on the BD, the white zone is the region for which the planet 
approaches the BD without falling on it. This last situation hap- 
pens for planets inside the corotation radius which are saved by 
the acceleration of the BD, however the BD tide is not strong 
enough to push them back farther than their initial position. The 
zone between the corotation line and the white zone corresponds 
to planets that begin falling but are then pushed farther than their 
initial position. 



Delta o/oq 
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Fig. 7. Map of Aa/a for a range of initial semi-major axis and different 
BD masses. Each BD has the same initial rotation velocity. The dashed 
dotted line is the corotation distance at the beginning of integration and 
the dashed line is the Roche limit. 



Figure 7 shows a map of Aa/ao as a function of the initial 
semimajor axis and BD mass. 

In contrast with Fig. 6 of the previous subsection, this map 
shows that the planets inside the corotation distance are effi- 
ciently pushed away. The initial corotation distance scales here 
as (Af* + M p ) 1/3 ~ M* 1/3 as we got rid of the initial BD rotation 
dependance. 

Within this set of calculations with mass-independent initial 
spin rates, it is interesting to note that the planets that evolve 
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according to the intermediate case - planets that start interior 
to corotation but survive - are more common around lower mass 
BDs. This is because of the combination of the small masses and 
small radii of low mass BDs: planets inside the corotation radius 
fall slowly and thus have a longer time window for the corotation 
radius to catch up with them. This tendency is absent in Fig. 6 
because the different initial rotation velocities of the BD cancels 
this effect. 

Furthermore, with identical initial spin rates low- and high- 
mass BDs push the planets farther than medium mass BDs. It 
is logical that higher mass BDs being heavier and bigger push 
planets farther but for low mass BDs the explanation is not as 
straightforward. The reason for this is that the low mass BDs 
evolve less than high mass BDs and they spin up relatively 
slowly because their radii do not decrease much. Surviving plan- 
ets remain close to the BD for a long time so the tides remain 
active and keep on pushing the planets outward on timescales of 
10 Gyrs. 

3.2.4. Influence of the time zero 

Gaseous protoplanetary disks are observed to dissipate after a 
few Myr around young stars (Haisch et al., 2001; Hillenbrand 
et al., 2008). Disks around low-mass stars have systematically 
longer lifetimes, of up to roughly 10 Myr (Pascucci et al., 2009). 
The evolution of discs around BDs is not well constrained but 
it is reasonable to assume lifetimes between one and ten mil- 
lion years. If planets form at the same approximate tempera- 
ture around BDs as around main sequence stars then their for- 
mation times should be shorter by a factor of 40-1000, since 
the accretion time scales inversely with the orbital frequency 
(see Safronov (1969); for a comparison of simulated and calcu- 
lated accretion times around low-mass stars see Raymond et al. 
(2007)). Thus, planets around BDs probably form before the disk 
dissipates and so we consider "time zero" for our tidal calcula- 
tions to be the time of the dissipation of the gaseous protoplane- 
tary disk. 

We now compare the tidal evolution of planets with time zero 
of 1 and 10 Myr, which should bracket the range of reasonable 
values based on the arguments presented above. The only impor- 
tant difference between these two initial conditions is between 
one and ten million years the BD contracts and spins up consid- 
erably. Thus, at 10 Myr the corotation radius is smaller (see Fig. 
1). 

We performed two suites of tidal calculations of 1M planets 
evolving around BDs with masses from 0.01 - 0.08M© for time 
zero values of 1 and 10 Myr. 

Maps of Aa/ao as a function of the initial semimajor axis and 
BD mass are presented in Figures 6 (for a time zero of 1 Myr) 
and 8 (for a time zero of 10 Myr). 

The corotation distance varies as £l~ 2 p 3 (M BD + M p ) 1/3 ~ 

a ~BD M BD SUCh that f0r 

our chosen initial BD rotation rates 
(see Fig. 1) the corotation distance increases with increasing BD 
mass. Consequently, the range of semi-major axis resulting on 
surviving planets decreases with mass. High-mass BDs cannot 
have planets as close-in as low-mass BDs. We can thus infer that 
there is a forbidden region larger than the one delimited by the 
Roche limit where no planet of 1M© should survive. This rough 
limit may be useful for future observations. 

Figures 6 and 8 also show that planets orbiting low-mass 
BDs are pushed farther than those orbiting high mass BDs. This 
is due to the smaller initial rotation period of low-mass BDs, 
they compensate for their smaller masses by being fast rotators. 



From these maps, we can also observe that Aa/ao decreases to 
as the initial orbital distance increases, which is logical because 
tidal effects decrease quickly with distance. 



Delta a/a 




0.01 0.02 0.03 0.04 0.05 0.06 0.070.08 
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Fig. 8. Map of Aa/a for a range of initial semi-major axis and different 
BD masses. Time zero is 10 7 yrs. The purple zone has shrunk compared 
to Fig. 6 because the BD has spun up and the corotation radius is thus 
smaller. 



After 1 Myr the BD's rotation period is 36 hours and by 
10 Myr it decreases to about 10 hours. During the interval from 
1-10 Myr the corotation radius has shrunk correspondingly so 
there is a wider range of initial semi-major axes that allows the 
planets to survive tidal evolution, which may be advantageous 
for planet detection (see Discussion in Section 4). 

The corotation distance does not have quite the same shape 
at 1 and 10 Myr (see Figs 6 and 8) because of the dependance 
on £Ibd- Indeed, we have seen that the time evolution of the spin 
rate varies with the BD mass (Section 2.3.3). Fig. 1 shows that 
the BD rotation period no longer increases monotonically with 
BD mass after 10 Myr, especially for BD masses larger than 
0.03 M Q . This explains the decrease in the corotation distance 
between M BD = 0.03 M© and 0.04 M© in Fig. 8. For BDs with 
masses below 0.16M©, the Roche limit is outside the corotation 
distance so some planets exterior to corotation do not survive 
the evolution. As shown in 8, the forbidden region is defined as 
the maximum of the Roche limit and the corotation distance. In 
contrast with a time zero of 1 Myr (Fig 6), here the corotation 
distance is coincident with the outer limit of the crash region. 

For to = 10 7 yrs the BD rotates faster than at to = 10 6 yrs 
so planets that survive are pushed away more strongly and the 
white zone where planets approach the BD without falling on it 
is even narrower. 

3.2.5. Influence of planetary mass 

Tidal effects are stronger for more massive bodies, so if the plan- 
etary mass is increased the planets are pushed farther or fall 
quicker. Moreover, the corotation radius increases as (M BD + 
M p ) 1/3 so heavier planets have slightly higher a SW i tc h. 

Figure 9 shows a set of tidal evolution calculations for plan- 
ets with masses from 0.1-10 M©, starting from just exterior to 
the corotation radius. For a planet of 0.1 M© tides only have the 
effect of reducing the eccentricity, but there is no semi-major 
axis change. The planet mass has little effect on the eccentric- 
ity evolution; in all cases e is damped to zero on a 1 - 10 Myr 
timescale. However, the final semi-major axis depends strongly 
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Fig. 9. Tidal evolution of planets with different masses starting at 9 x 
10~ 4 AU around a 0.04M© BD. The heavier planets are pushed away 
farther. The planetary masses are : 0.1, 0.5, 1, 2, 5, 10 M e . 



on the planet mass. This is because, especially at early times be- 
fore the BD contracts, the BD tide is stronger for larger planet 
masses and this is what drives the orbital expansion. 

For completeness, we ran a limited set of tidal calculations 
with planet masses of one Saturn mass and one Jupiter mass. 

The planet's dissipation factor was taken to be the same one 
as the BD - so the same as Hansen (2010)'s cr p . What is inter- 
esting for the case of gas giant planets is that they can have a 
significant influence on the tidal evolution of the BD. When the 
massive planet is close, angular momentum is transferred from 
the BD's rotation to the planet's orbit much more efficiently than 
for terrestrial planets. 

Figures 10 and 1 1 show the evolution of the gas giant's semi- 
major axis, eccentricity, and rotation period and the ratio of BD 
to planetary tides for a Saturn- and Jupiter-mass planet. The 
semimajor axis evolution for gas giants is qualitatively similar 
to the evolution for lower-mass planets, with larger- amplitude 
outward migration due to the planets' larger masses. But this is 
where the similarities end. 

The eccentricity evolution is dramatically different for gas 
giants: the planetary tide still acts to decrease the eccentricity 
but the BD tide can act to increase it. The phase diagram in Fig. 
1. of Leconte et al. (2010) (reproduced in Fig. 12) shows that 
for high Qbd/k the BD tide can indeed increase the eccentricity. 
This can be interpreted as a slingshot effect: the planet is kicked 
when approaching periastron, and if the kick is strong enough 
the trajectory gets more and more eccentric (Hut, 1981). 

At early ages the planets closest to the BD have a high mean 
motion and the BD has not accelerated yet so CIbdIk < jf ^(ej • 
This criterion comes from the BD tide contribution in Eqn. 3 and 
defines the limit e = in Fig. 12 (reproduced from Leconte et al. 



(2010)). For £l B D/n < 



18 Nel(e) 



the BD tide acts to decrease the 



1 1 Ne2(e) ' 

eccentricity. So for close planets at early ages, both BD and plan- 
etary tides act to decrease the eccentricity. This effect is visible 
in Figures 10 and 1 1 for ages between 10 3 yrs and 10 5 yrs. When 
the BD spins up, the criterion £l BD /n < jy ^(g) * s no l° n S er met 
so the BD tide acts to increase the eccentricity - this is the sling- 
shot effect. The planetary tide is not strong enough to counteract 
the eccentricity pumping from the BD. Figures 10 and 11 show 
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Fig. 10. Tidal evolution of a Saturn mass planet starting at different ini- 
tial semi-major axis around a 0.04 M© BD. The first panel shows the 
evolution of the semi-major axis (bold solid lines). The dashed dotted 
line represent the radius of the BD, the dotted line the corotation radius 
and the bold dashed line the Roche limit. The habitable zone is also 
plotted. The second panel shows the evolution of the eccentricity. The 
third panel shows the evolution of the rotation period of planets (solid 
lines) and the BD (dashed lines). The solid black line represents the ro- 
tation period of the BD if its evolution is governed only by conservation 
of angular momentum. The fourth panel shows the ratio of stellar over 
planetary tides and the horizontal line is a ratio of 1. 
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Fig. 11. Tidal evolution of a Jupiter mass planet starting at different ini- 
tial semi-major axis around a 0.04 M© BD. The third panel shows that 
a close-in Jupiter mass planet significantly slows the BD with respect 
to the rotation rate it would have if there was conservation of angular 
momentum. 



this phase between the ages of 10 5 yrs and 10 7 yrs. As the BD 
radius decreases, the BD tide weakens and after a few Gyr, the 
planetary tide becomes dominant and acts to once again decrease 
the eccentricity. 

Planets that start farther from the BD have slower mean mo- 
tion so SIbdIk > |f ^(g) anc * m e BD tide acts to increase the 
eccentricity. These planets skip the first stage of the evolution 
previously described. Their eccentricity is pumped by the BD 
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Fig. 12. cj/n vs e phase space, reproduced from Leconte et al. (2010) 
Fig. 1. Here, co = Q p or Q B d- The planets are in pseudo- synchronization 
so Q p /n remains on the black solid line, thus the planetary tide always 
contribute in making the semi-major axis and eccentricity decrease. If 
the eccentricity is ~ 0, we can see that if O^BD/n > 1 the BD tide con- 
tributes in pushing planets away and if Qbd/k > 18/11 ~ 1.64, the 
eccentricity increases. 



The key difference in evolution for gas giants come from the 
ability of massive planets to affect the BD's rotation such that the 
BD rotation period can no longer be approximated by the simple 
conservation of angular momentum during the BD's contraction. 
Gas giants are so massive that they act to reduce the acceleration 
of the BD significantly during contraction. Planets outside the 
corotation radius have slower mean motion than the BD's rota- 
tion rate, so the system tends to transfer angular momentum from 
the BD to the planet making the BD slow down and the planet 
spin up and therefore move outward. The stellar tide efficiently 
pushes the planets away. The efficiency of the spinning down 
of the BD is measured with f ti( i es (see Equation 11). This pa- 
rameter is negligible for terrestrial planets but for a Jupiter mass 
planet T BD is smaller so f tides is larger. Fig. 11 shows that for 
close-in Jupiter mass planets the rotation rate of the BD does not 
quite follow the evolution governed by the conservation of an- 
gular momentum. As in the first few million years the BD does 
not contract much, the effect is visible, but at t ~ 10 7 yrs the 
contraction of the BD becomes significant and the effect is over- 
powered and the rotation rate of the BD increases accordingly to 
the conservation of angular momentum. 



and when the BD tide gets weaker than the planetary tide the ec- 
centricity decreases. For planets farther than 0.04 AU, both tides 
are negligible and there is no effect on the semi-major axis or 
eccentricity. 

Figure 13 shows the evolutionary tracks of Jupiter mass plan- 
ets beginning at different semi-major axis in the phase diagram 
£l/n vs eccentricity. It shows the different phases of the tidal 
evolution as explained in the previous paragraph depending if 
£l BD /n is bigger or smaller than jf 

In this case, as can be seen in the fourth panel of Fig. 11, 
the BD tide dominates the evolution. So Fig. 13 shows only the 
effect of the BD tide on the evolution of the system except at the 
very end of the evolution when the planetary tide is becoming 
dominant. This is why there is a decrease in eccentricity in the 
e > region at the end of the evolutionary tracks. 
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Fig. 13. Evolutionary tracks of Jupiter mass planets orbiting a 0.04 M Q 
BD. The evolution starts on the solid black line of e = 0.01. The differ- 
ent solid lines are for different initial semi-major axis. The dashed line 
is where a = and the dashed dotted line is where e = 0. The quantities 
a and e are positive above the zero value line. 



3.2.6. Influence of BD dissipation 

The BD dissipation is a poorly-constrained parameter. Here we 
adopt Hansen's (2010) value for BDs and assume that it remains 
constant in time and constant for all BD masses. We note that 
Hansen (2010) considers that the dissipation factor for gas gi- 
ants is approximatively the same as for BDs. Indeed, gas giants 
and BDs are thought to have similar internal structures Baraffe 
et al. (2003b). As the BD can have a low rotation period it is 
interesting to compare the value of the dissipation factor we 
chose in comparison to a low rotation period gas giant dissi- 
pation factor. Leconte et al. (2010) gives a estimation of the 
dissipation in Jupiter k^jt^Tj < 2 - 3 x 10" 2 s, which corre- 
sponds to ctj ~ 6 x 10" 52 kg _1 m" 2 s _1 . Our dissipation factor 
being cr BD = 2.006 x 10" 53 kg _1 m~ 2 s _1 , the parameter space ex- 
ploration includes BDs which have the same dissipation factor 
as Jupiter. 

To test the effect of the BD dissipation on the tidal evolution, 
we performed a set of tidal calculations varying the dissipation 
factor cr BD by several orders of magnitude. In each integration, 
the BD mass was 0.04 M Q , the planet mass was 1 M©, and the 
planet started just exterior to the corotation radius. 

Figure 14 shows the results of tidal calculations varying cr BD 
from 0.001 to 1000 times our fiducial value. 

The higher the BD dissipation factor, the stronger the tides 
created by the planet in the BD, so the farther the planet is 
pushed. For BD dissipation factors of 0.05 x <tbd or smaller, 
tides are not strong enough to significantly change the planet's 
orbit. The eccentricity evolution depends on the BD dissipation 
factor. For a low dissipation factor, the eccentricity evolution is 
governed by the planetary tide and decreases to zero. For a high 
BD dissipation factor - for dissipation factor > 10 cr BD - a signif- 
icant change can be seen in the eccentricity evolution as shown 
in Figure 14. The effect is even more significant for a BD dis- 
sipation of 1000 cr BD , indeed the eccentricity decreases quickly 
during the first 10 4 yrs because both tides contribute to decreas- 
ing the eccentricity. At t = 10 4 yrs, £l B o/n becomes larger than 
if Afe2(g) ' wmc h defines the e = line in a phase diagram. So the 
eccentricity increases and this effect is visible because the BD 
dissipation factor is high enough to make the BD tide dominant. 
At t ~ 10 7 yrs, the BD starts spinning up significantly, so the 
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BD tide becomes negligible and the planetary tide contributes in 
making the eccentricity decrease once more. 
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Fig. 14. Tidal evolution of a 1 M e planet starting at 9.5 x 10~ 4 AU around 
a 0.04 M Q BD for different BD dissipation. The planet is pushed farther 
if the BD dissipation is higher. The BD dissipation factors are : 0.001, 
0.05, 0.1, 1, 10, 20, 1000xo- 5D . 



3.2.7. Influence of the planetary dissipation 

We calculate the planetary dissipation using the estimate of 
Neron de Surgy & Laskar (1997) for the tidal lag time ki&At®. 
Eqn (6) allows us to obtain the planetary dissipation factor as 
used by (Hansen, 2010). The planetary dissipation factor also 
depends strongly on the planetary radius: for a given mass and 
tidal dissipation factor, a larger planet experiences a stronger 
tidal flexion and dissipates more heat. 

The planetary tide is important when the eccentricity is non 
negligible. As the planet is in pseudo-synchronization it is con- 
stantly under changing stress which brings about more tidal dis- 
sipation. The system evolves to decrease the eccentricity to min- 
imize the dissipation. 

Figure 15 shows the influence of the planetary dissipation 
on the tidal evolution of a 1 M® planet orbiting a 0.04 M BD. 
The eccentricity decreases faster when the planetary dissipation 
factor is higher because the tides act more quickly to minimize 
the stress in the planet. Indeed, when planetary dissipation in- 
creases, the dissipation timescale T p decreases. Within just 1000 
years significant differences can be seen between planets with 
different values of cr p . From an initial eccentricity of 0.01, the 
planet with the largest dissipation factor had its eccentricity cut 
in half in 1000 yrs, while the planet with the lowest dissipation 
factor needs ~ 10 8 yrs for this to occur. 

Planets with higher planetary dissipation rates are pushed 
farther from the BD. This is due to the fact that for low plan- 
etary dissipation, the eccentricity remains non negligible for a 
long time. Because of this eccentricity, the planetary tide coun- 
teracts a portion of the BD tide and thus the planet is not pushed 
as far. However, this is a relatively weak effect because changing 
the planetary dissipation by several orders of magnitude barely 
affects the final orbital distance. Indeed, for a planet beginning 
its evolution at an initial semi-major axis of 9x 10" 3 AU, the final 
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Fig. 15. Tidal evolution of a 1 M e planet starting at 9 x 10" 4 AU around 
a 0.04 M Q BD for different planetary dissipation. Dependance on the 
planetary dissipation. The planet is pushed farther if the planetary dissi- 
pation is higher. The planetary dissipation factors are : 0.01, 0.1, 1, 10, 
100 xcr D 



positions differ by only 2.5% after 10 Gyrs. If the planet initial 
semi-major axis is 12 x 10" 3 AU, they differ by only 0.5%. 

The planetary dissipation factor has little influence on the fi- 
nal semi-major axis for low initial eccentricity. The fact that the 
planet's orbit expands is almost entirely due to the BD tides. This 
has important implications for observations because if a planet 
is discovered orbiting a BD and its mass and orbit are known, 
then we can infer a value for the BD dissipation. However the 
difficulty in making that assessment is that it is only true for a 
system in which the eccentricity remained small throughout the 
evolution. Observing at time t does not tell us about the evolution 
of the eccentricity. A combination of the initial eccentricity and 
the planetary dissipation can be very roughly estimated in sys- 
tems with multiple planets (e.g., (Papaloizou & Terquem, 2010; 
Barnes et al., 2009a)) but that is beyond the scope of this paper. 

Figures 17 and 16 show the effect of varying both the initial 
eccentricity and the BD (Fig. 17) and planetary (Fig. 16) dissi- 
pation factors. The initial eccentricity has little influence on the 
dependance on BD dissipation because the eccentricity evolu- 
tion is mainly determined by planetary tides. The eccentricity 
decreases to near zero in 10 6 yrs to a few times 10 7 yrs. 

However, varying the initial eccentricity has a strong influ- 
ence on the effect on planetary dissipation. The eccentricity evo- 
lution itself is straightforward: e decreases to negligible values 
in a few 10 4 yrs for the highest planetary dissipation factors to 
10 9 yrs for the lowest dissipation factors. Higher eccentricity de- 
creases faster than smaller ones so that for a same planetary dis- 
sipation factor both eccentricities reaches low values at about the 
same time. 

What is particularly interesting is that the initial eccentric- 
ity influences the final semi-major axis for different planetary 
dissipation rates. For low planetary dissipation, the eccentricity 
remains non negligible for a long time and the evolution is com- 
plicated. The planetary tide is weak so the BD tide acts to slowly 
push the planet away, but when after a few million years the 
BD has shrunk significantly the BD tide's importance decreases 
to the point where the planetary tide eventually wins over and 
pushes the planet inwards. 
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Fig. 16. Tidal evolution of a 1 M e planet starting at 9 x 10~ 4 AU around 
a 0.04 M Q BD for different BD dissipation and eccentricity. The BD dis- 
sipation factors are : 0.05, 0.1, 1, 10, 20 x <x* and the initial eccentricity 
0.01 and 0.2. 
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Fig. 17. Tidal evolution of a 1 M e planet starting at 9 x 10~ 4 AU around 
a 0.04 M Q BD for different planetary dissipation and eccentricity. The 
planetary dissipation factors are : 0.01, 0.1, 1, 10, 100 x cr p and the 
initial eccentricity 0.01 and 0.2. 



If the planet is in pseudo-synchronization, the contribution 
of the planetary tide always pushes the planet inward and cir- 
cularizes the orbit (see Fig. 12). The only difference between 
low dissipation and high dissipation planets is then the timescale 
of the evolution. For higher planetary dissipation, the planetary 
tide is strong enough to quickly decrease the eccentricity and 
the evolution is then determined by the BD tide, whereas for 
lower planetary dissipation the eccentricity remains high for a 
long time and only at late times, when the BD tide becomes neg- 
ligible, is the evolution determined by the planetary tide. 



4. Discussion 

4.1. Brown dwarfs and habitable zone 

We now consider the habitability of planets orbiting BDs in 
the context of the circumstellar habitable zone, i.e. the range 
of orbital radii for which a planet's surface could have liquid 
water (Kasting et al., 1993; Selsis et al., 2007). The issue of 
habitability is particularly interesting in the context of BDs be- 
cause, as the BD contracts and cools the habitable zone moves 
inward. Meanwhile, in most cases tidal evolution push plan- 
ets outward. Given this divergent radial evolution, planets that 
survive around BDs pass through the habitable zone and will 
thus spend a finite time with potentially habitable conditions. 
Andreeshchev & Scalo (2004) found that for BD masses larger 
than 0.04 M© a planet can remain in the habitable zone for more 
than 4 Gyrs if it forms sufficiently close to the Roche limit. 
However, Andreeshchev & Scalo (2004) did not take tidal evo- 
lution into account. As we will see below, tidally-driven outward 
migration reduces this habitable time considerably. 

We calculated the location of the habitable zone in two dif- 
ferent ways. We first assumed that a planet on the outer edge 
receives 400 W/m 2 , the same flux as the ~ 8M e potentially- 
habitable exoplanet GJ 581d (Udry et al, 2007; Wordsworth 
et al., 201 1). The inner edge of the habitable zone is taken to be 
four times this flux, 1600 W/m 2 ; although this is not very phys- 
ical it offers a reasonable approximation. The 400 - 1600 W/m 2 
habitable zone is marked with solid lines in Figs. 2, 9, 10, 11, 14, 
15,16 and 17. Here the most important limit is the outer limit be- 
cause it determines the time a planet spends on the HZ. For our 
second estimate we used the same method but assumed that a 
planet on the outer edge receives 240 W/m 2 , which is thought 
to be the minimum flux a planet requires to be habitable with 
the best atmospheric conditions possible. This second habitable 
zone estimate, of 240 - 480 W/m 2 - is shown with the dashed 
line in Figs. 2, 9, 10, 11, 14, 15, 16 and 17. We note that, given 
the diverse ways to heat a planet (e.g., Forget & Pierrehumbert 
(1997), Pierrehumbert & Gaidos (2011)), the outer edge of the 
habitable zone is in reality a function of at least planet's atmo- 
spheric properties. 

We used the same tidal calculation presented in section 3 to 
study the influence of different parameters on how much time a 
planet spends in the habitable zone (see Figs. 9, 14 and 15). 

Planets at larger orbital distances spend less time in the hab- 
itable zone because it sweeps past them early in their evolu- 
tion. Increasing the planet mass reduces the time spent in the 
habitable zone (Fig. 9 ). Lighter planets are less influenced by 
tides so they will tend to have an almost constant semi-major 
axis. For a given initial semi-major axis, lighter planets spend 
more time in the habitable zone than heavier planets. Moreover, 
if planets form with substantial eccentricity they receive more 
flux early in their evolution, which is not in favor of habitabil- 
ity. Planets in their early evolution are inside the inner edge of 
the HZ so they are already too hot for liquid water, a non zero 
eccentricity will just worsen the situation. For standard dissi- 
pation factors and for planets initially close to their BD, the 
eccentricity drops to near-zero before the planet reaches the 
habitable zone, or rather before the habitable zone reaches the 
planet. However, there do exist cases where the planet starts far- 
ther from the BD so the eccentricity remains high for a long 
time, which can contribute to keeping the planet habitable for 
a slightly longer time. Indeed, the flux received by the BD is 
F(a, e ^ 0) = F(a, e = 0)(1 - e 2 )~ 112 (see also Barnes et al., 
2009b). However, the flux increase is only significant for very 
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eccentric planets (e > 0.8), and such planets would need to be 
relatively distant from their BD to maintain a high eccentricity 
high for a long time. 

Increasing the BD dissipation and the planetary dissipation 
also decreases the time a planet spends in the habitable zone be- 
cause the planets are pushed farther and the HZ sweeps past them 
early in their evolution. However, the effect of the BD dissipa- 
tion is more pronounced than of the planetary dissipation, which 
is seen in Figs. 14 and 15. The time spent in the habitable zone 
decreases faster if we increase the BD dissipation factor rather 
than the planetary dissipation factor. 

The main parameter that determines a planet's habitable zone 
lifetime is the mass of the BD. More massive BDs are more 
likely to host a habitable planet than low mass BDs. This is due 
to the evolution of the luminosity as a function of BD mass. The 
luminosity of high mass BDs is higher and decreases less than 
for low mass BDs. So the HZ for high mass BDs is farther from 
the BD and this is favorable for habitability because tidal evo- 
lution tends to be weaker such that outward migration is slower 
and less pronounced. 




0.01 0.02 0.03 0.04 0.05 0.06 0.07 0.08 
BD mass (M Sun ) 



Fig. 18. Time spent in the HZ by a 1 M e orbiting a 0.04 M Q BD. Time 
zero is 1 Myr. This time has been calculated with the most optimistic 
case, which is for a switch . 



To sum up, the most optimistic scenario for habitability ex- 
ists for BDs with masses above 0.04 M© and for parameters that 
induce the smallest orbital change, that is low BD and planetary 
dissipation and low planetary mass. However, it is important to 
keep in mind that the final position of planets and thus the time 
they spend in the HZ depends predominantly on the initial semi- 
major axis compared to the corotation radius. Figure 19 shows 
that for low BD dissipation, the best case corresponding to a SW i tc h 
leads to orbital changes and actually brings the planet closer to 
the BD than it was initially. 

4.2. Observational possibilities 

BDs are thought to be numerous and their small masses and 
radius makes them better targets for radial velocity and transit 
searches than main sequence stars for a given limiting magni- 
tude (Blake et al, 2008). 

Observing planets orbiting BDs could help constrain tidal 
dissipation models. In Section 3 we showed that the final or- 
bital distance of a planet depends mainly on its initial position 
compared to the initial corotation radius. It is therefore hard to 



predict anything from a single observation of a planet orbiting 
a BD. However, a large sample of planets orbiting BDs could 
place constraints on certain tidal parameters. 



100 r 




10 3 10 4 10 5 10 6 10 7 10 8 10 9 10 1 



time (yrs) 

Fig. 19. Tidal evolution of the rotation period a 1 M e planet orbiting 
a 0.04 M Q BD. The different panels are for different BD dissipation 
factors : 0.05, 0.1, 1, 10, 20 X cr BD . The initial eccentricity is 0.01. 



The main quantity that could be constrained with a sample 
of planets orbiting BDs is the ctbd, the dissipation rate within 
the BD. Figure 19 shows the evolution of the rotation period of 
a suite of 1 M® planets orbiting a 0.04 M BD. The eccentricity 
drops to zero and the planet is in true spin-orbit synchronization, 
so the y-axis in Fig. 19 represents the orbital period. As we saw 
in subsection 3.2.6, the BD dissipation factor strongly influences 
the final planetary orbit. 

Fig. 19 shows that for low BD dissipation the distribution 
of final orbital distances is less densely packed than the distri- 
bution of initial orbital distances, with orbital periods as small 
as ~ 10 firs. The planets' evolution is dominated by planetary 
tides, which bring the planets closer to the BD and effectively 
spread the initial planet distribution. In contrast, for the highest 
BD dissipation the distribution of final orbital distances is much 
more densely packed than the distribution of initial orbital dis- 
tances and at much larger orbital periods because the BD tide 
dominated and planets have all been pushed outward. The plan- 
ets accumulate at an orbital period of ~ 80 hrs. 

Figure 20 shows the minimum final orbital period for 1 M e 
planets as a function of the BD mass and the BD dissipation. As 
expected, the minimum period increases for higher BD dissipa- 
tion. For low BD dissipation, the minimal period is about 20 hrs 
and does not change much with mass. The fluctuations observed 
are due to the extreme sensitivity to initial conditions because the 
final orbital distance of the planets depends on when they cross 
the shrinking corotation distance. However, the decrease of the 
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Fig. 20. Final minimal period of 1 M e planets as a function of the BD 
mass and for different BD dissipations. 



final rotation period from a BD mass of 0.01 M Q to 0.02 M© 
was explained in subsection 3.2.2 and is due to the fact a BD of 
0.01 M Q pushes the planets during its entire evolution. For dissi- 
pations of 1, 10 and 20 ctbd, the minimal final period increases 
between masses of 0.02 M© and 0.07 M© because the more mas- 
sive the BD, the stronger the BD tide. However the final decrease 
between 0.07 M© and 0.08 M© is due to sensitivity to initial con- 
ditions. For example, for M* = 0.07 M©, the planet crosses the 
corotation radius at a few 10 5 yrs and for M* = 0.072 M© it 
crosses it at ~ 10 6 yrs. So the planet around the 0.07 M© BD is 
pushed outward sooner when the BD radius is still large so the 
planet is pushed more strongly and for a longer time. The preci- 
sion to which a switch has been determined here is not sufficient to 
remove these fluctuations . 

Thus, with a sufficiently large sample of planets orbiting BDs 
with well-characterized masses and orbits, it should be possible 
to infer - or at least constrain - values for the dissipation in BDs. 

Figures 19 and 20 assumed a time zero of 1 Myr. However, 
as seen in Section 3.2.4, tidal evolution is weaker if time zero 
is later simply because BDs contract significantly during 1 and 
10 Myr. Planets can form and survive closer to their BD hosts if 
time zero is later, i.e. for longer lifetimes of the gaseous proto- 
planetary disk. Thus, planets may have even shorter orbital peri- 
ods than those shown in Fig. 20 for later time zero. For a 1 M® 
planet orbiting a 0.04 M© BD and a time zero of 10 Myr, the 
final minimal periods are of about 10 hr for BD dissipation of 
0.05 and 0.1 cr BD , 20 hr for 1 cr BD , 30 hr for 10 cr BD and 40 hr 
for 20 <tbd- If time zero is 10 Myr, the final minimal period is 
half it is if time zero is 1 Myr. 



5. Conclusion 

In this paper, we calculated the tidal evolution of planets around 
BDs using a standard constant time lag tidal model (e.g., Leconte 
et al 2010, Hansen 2010). Our calculations account for the fact 
that BDs contract with time (Baraffe et al., 2003b) and therefore 
spin up. We calibrated the BD initial rotation period with ob- 
servational data (Herbst et al., 2007). We then investigated the 
parameter space in order to gauge each parameter's effect on the 



evolution of planets orbiting BDs and on the time planets can 
spend in the habitable zone. 

We found that the most important parameter was the initial 
orbital distance compared to the corotation distance. A planet 
outside the corotation radius migrates outward and stops at a dis- 
tance where tidal effects are negligible. A planet inside the coro- 
tation radius evolves inward and is likely to eventually merge 
with the BD. However, as the BD spins up with time the corota- 
tion radius shrinks and in some cases the corotation radius can 
pass by the planet orbit and reverse the direction of the migra- 
tion. The fact that BDs contract with time allows the system to 
experience many different evolutions given the importance of pa- 
rameters such as the initial eccentricity and the dissipation fac- 
tors of the planet and the BD. The planetary tide and the BD tide 
vary in intensity and the value of their ratio is extremely helpful 
to understand which behavior is due to which tide. 

Thus the study of tidal evolution of a planet orbiting a BD 
is interesting from a theoretical point of view because of the va- 
riety of possible evolutions due to the spin-up of the BD and 
the effect of parameters like the dissipation factors. But it is also 
interesting from an observational point of view because planets 
that form close to their BD can have very short final orbital pe- 
riods. In the most favorable cases, their periods can be as short 
as ~ 10 hrs. Those planets could be observed in transit around 
BDs. 

For BDs more massive than 0.05 M , planets can spend a 
significant amount of time in the habitable zone - from 1 Gyr to 
-10 Gyrs- so this raises also the question of habitability of such 
planets. Furthermore, planets with pseudo- synchronous or syn- 
chronous rotation are interesting case studies for climate model- 
ing (Wordsworth et al., 2011). 
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